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RESUME. Une modelisation 3D et multi-echelles du comportement physico-chimique et 
mecanique du materiau cimentaire est developpee. L’heterogeneite du materiau est 
representee dans des volumes elementaires representatifs des differentes echelles d’ etudes. 
Les proprietes physico-chimiques et mecaniques a une echelle se determinent a partir de 
cedes de V echelle inferieure. U echelle de la micro structure, id V echelle des hydrates est 
modelisee en utilisant le code CEMHYD3D. Ce logiciel est ensuite couple au logiciel 
ABAQUS pour le calcul des proprietes mecaniques. La demarche est ainsi appliquee pour 
V etude du comportement mecanique d’un mortier a E/C=0, avant et apres lixiviation par 
Veau. Les resultats des simulations sont en bon accord avec ceux de V experimentation. 

ABSTRACT. 3D multi-scale and multi-physics modeling of cement based materials is developed. 
Two numerical tools are combined to predict the mechanical behavior of these materials 
CEMHYD3D and ABAQUS. The outcomes of the simulation at the micro-scale are used at 
the meso-level modeling. This approach is then applied to investigate the mechanical 
behaviour of a mortar with water-to-cement ration equal to 0.4. The effect of the leaching 
phenomenon is studied. The numerical results of the modeling are consistent with the 
experimental ones. 

MOTS-CLES : modelisation multi-echelles, micro structure, comportement mecanique, 

lixiviation, ciment. 
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Introduction et contexte 

La durabilite des ouvrages en beton est une preoccupation de plus en plus mise 
en avant par les gestionnaires d’ ouvrages. Ces derniers exigent des durees de vie de 
plus en plus longues avec un minimum d’entretien. Or la durabilite des ouvrages en 
beton est etroitement liee a celle du beton lui-meme. Ce dernier est un materiau 
heterogene et poreux en equilibre chimique reagissant avec son environnement. 
L’ etude de sa durabilite fait intervenir des couplages entre plusieurs disciplines : 
chimie, physique et mecanique, ceci a diverses echelles de la nanostructure a la 
macrostructure. Les etudes experimentales des phenomenes lies a la durabilite 
necessitent souvent un temps tres long et les resultats sont parfois difficiles a 
exploiter. La simulation numerique est un outil precieux qui permet un gain de temps 
consequent mais aussi une meilleure comprehension des phenomenes mis en jeu. La 
prediction du comportement physico-chimique et mecanique du beton et de la 
structure au cours du temps necessite le developpement et 1’ utilisation de modeles 
couples. La robustesse de la modelisation est etroitement liee a la qualite de la 
representativite du materiau numerique lui-meme et des phenomenes mis en jeu. 

Les travaux menes ici visent a developper une modelisation tridimensionnelle et 
multi-echelles du comportement physico-chimique et mecanique du materiau 
cimentaire. La modelisation multi-echelles et multi-physiques est rendue aujourd’hui 
possible grace aux progres informatiques et au developpement d’outils numeriques 
puissants. L’approche developpee consiste a representer I’heterogeneite du materiau 
dans des volumes elementaires representatifs des differentes echelles d’ etudes 
(micro-meso). Les proprietes physico-chimiques et mecaniques a une echelle se 
determinent a partir de cedes de 1’ echelle inferieure. L’ echelle de la microstructure, 
ici r echelle des hydrates (quelques microns) est modelisee en utilisant le code 
CEMHYD3D developpe par le NIST/USA (Bentz, 1997). Ce logiciel est ensuite 
couple a un autre logiciel, ABAQUS, dedie aux calculs mecaniques et de transferts. 

Nous presenterons d’abord la demarche et la methodologie adoptee. Ensuite, 
nous I’appliquerons pour I’etude d’un mortier a E/C = 0,4. Nous nous interesserons 
plus particulierement au comportement mecanique du materiau sous traction et 
compression et ensuite a sa durabilite vis-a-vis de la lixiviation par I’eau pure. Les 
resultats des simulations seront confrontes aux resultats experimentaux. 

2. Descriptif de la demarche et de la methodologie 

La demarche adoptee ici comporte deux etapes principals, schematisee dans la 
figure 1 et decrite ci-dessous : 

Etape 1 : Representation numerique en 3D du materiau etudie 

Le beton est un materiau composite heterogene compose d’un squelette 
granulaire enrobe dans une matrice poreuse constitute par la pate de ciment 
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hydratee. La granulometrie du squelette granulaire varie du micron a la dizaine de 
millimetres. Trois echelles sont necessaires pour decrire son comportement : 

Echelle macroscopique : le materiau est modelise comme un empilement 
granulaire 3D compact de grains de dimensions comprises entre quelques mm et 
quelques dizaines de mm (gravillon), enrobes par un mortier. La zone d’ interface 
entre les grains et la matrice appelee aureole de transition a des caracteristiques 
mecaniques et des proprietes de transfert differentes de celle de la pate de ciment du 
mortier. Les proprietes du mortier sont deduites de la mesostructure generee a 
r echelle mesoscopique. 

Echelle mesoscopique : le mortier est decrit comme un empilement compact 3D 
de grains de dimensions comprises entre une centaine de microns et quelques 
millimetres (sable), enrobes par une pate de ciment durcie. Les proprietes de la pate 
de ciment sont deduites de la micro structure generee a 1’ echelle micro scopique. 

Echelle microscopique: la pate de ciment durcie est decrite a partir de la 
composition du ciment et du rapport E/C comme un assemblage poreux 3D 
compose: 

- de phases hydratees telles que les silicates de calcium hydrates, la portlandite, 
les sulfo-aluminates de calcium hydrates (ettringite, monosulfoaluminate), etc... 

- de phases anhydres ou en cours d’hydratation telles que les silicates de 
calcium et les aluminates de calcium. 


Etape 1 : Modelisation en 3D du materiau par CEMHYD3D 


Echelle micro-scopique Echelle meso-scopique Echelle macro- scopique 
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Etape 2 : Calcul des proprietes mecaniques et de transferts par 
ABAQUS 


Figuel. Schemas illustrant la demarche de moddisation multi-echelles adoptee 
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Etape 2 : Calcul des proprietes mecaniques et de transfer! du materiau sain ou 
altere 

Le materiau 3D genere dans 1’ etape precedente est ici repris par un autre logiciel 
plus dedie a la mecanique et a 1 ’ etude des phenomenes de transferts : ici le logiciel 
ABAQUS. L’image pixelisee est alors transformee en un maillage Elements Finis 
utilisable par ce logiciel. Chaque voxel representant une phase du materiau est alors 
considere comme une maille cubique dans ABAQUS. L’ alteration pent aussi etre 
directement prise en compte dans cette etape. 


3. Application de la methodologie a Fetude d’un mortier 

3.1. Materiau 

La methodologie developpee plus haut a ete appliquee a un mortier qui a fait 
I’objet d’une etude experimentale anterieure (Kamali, 1999) et (Le Bellego, 2001) 
afin de caracteriser son comportement sous I’effet de sollicitation de compression, 
traction, flexion avant et apres lixiviation par du nitrate d’ ammonium. Sa 
composition ainsi que celle du ciment le constituant sont decrites dans ces 
references. Ce mortier a base de ciment portland CEM I 42,5 a un rapport E/C egal a 
0,4. Notre etude sera done limitee aux echelles micro et meso. 

3.2. Etude du materiau sain 

3.2.1. Modelisation de la micro structure de la pate de ciment et de la mesostructure 
du mortier 

Un cube de pate de ciment portland a E/C=0,4 de 100 x 100 x 100 voxels est 
simule, chaque voxel correspond a 1 pm^. La figure 2 presente une coupe 2D de 
I’image 3D representant la pate hydratee. Les quantites d’hydrates et la valeur de la 
porosite capillaire sont donnees dans le tableau 1 . 


Tableau 1. Teneur en hydrates et porosite capillaire de la pate de ciment hydratee 
d’ apres CEMHYD3D 


Porosite 

CH 

C-S-H 

(C 3 AH 6 , AFm et AFt) 

17% 

14,7% 

40,2% 

17,9% 


La mesostructure du mortier est representee par un cube elementaire de 5 mm de 
cote contenant une quantite XI de spheres de diametre 1,125 mm, une quantite X2 
de spheres de diametre 0,125 mm et une quantite X3 de spheres de diametre 0.375. 
Les quantites, XI, X2 et X3 s’obtiennent a partir respectivement des teneurs 
volumiques en sable SI, S2 et en bulles d’air. Autour des gros grains de diametre 
1,125 mm, il est possible de prendre en compte la presence d’une aureole de 
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transition. La figure suivante presente une coupe 2D d’un cube de mm de cote issu 
de la mesostructure 3D du mortier simule. 


Composants 
I I du ciment 
I — I non hydrate 
Porosite 
□□ C-S-H 
H CH 
□□ AFt 

□□ C3AH6 




Figure 2. Images 2D de la micro structure 3D de la pate de ciment hydratee 
(gauche) et image 2D de la mesostructure 3D du mortier simule (droite) 


3.2.2. Calcul des proprietes mecaniques du mortier 

Pour des raisons de capacite de memoire, les calculs mecaniques avec ABAQUS 
ont ete realises sur des cubes elementaires de 50 qm de cote pour la pate de ciment 
et 5 mm pour le mortier. D’apres Kanit et al (2003), la representativite d’un volume 
elementaire peut-etre verifiee en imposant differentes conditions aux limites 
(conditions homogenes en deplacements ou en forces). Ce critere a ete verifie pour 
les volumes elementaires consideres dans cette etude. L’ aureole de transition n’est 
pas prise en compte. 

Les caracteristiques de la pate de ciment hydratee ont d’abord ete calculees. Une 
maille cubique represente les dimensions d’un voxel de Iqm^. A chaque phase de la 
pate de ciment ont ete attributes des caracteristiques mecaniques : module elastique, 
coefficient de Poisson, resistance a la traction. Les valeurs de module elastique et de 
coefficient de Poisson des differentes phases considerees sont obtenues par nano- 
indentation. Elies sont issues de la litterature et ont deja ete presentees auparavant 
(Kamali, 2004). La valeur de la resistance a la traction d’une phase est consideree 
proportionnelle au module elastique (ft = E/10000). Un critere de Rankine est 
considere pour la rupture des differentes phases. A I’echelle de la mesostructure 
(mortier), un modele de fissuration diffuse (Hillerborg, 1976) est utilise pour prendre 
en compte le comportement en traction de la pate de ciment hydrate. Les parametres 
de ce modele (resistance en traction, courbe contrainte-deplacement post-rupture) 
sont identifies par calage sur la courbe contraintes-deformation en traction du VER 
de pate de ciment. Les grains de sable suivent le meme type de comportement. 

Par ailleurs, a cette echelle, la rupture est supposee etre la consequence d’etats de 
contraintes de traction locaux generes par I’heterogeneite du materiau. Ainsi la seule 
connaissance du comportement en traction de la pate de ciment est suffisante pour 
obtenir a la fois la reponse du mortier a une sollicitation de traction et de 
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compression. La figure 3 montre les courbes contraintes-deformations du mortier et 
permet de conclure que les resultats de la modelisation sont proches des resultats 
experimentaux. Les facies de fissuration sont egalement bien representes (figure 4). 



Figure 3. Comportement sous chargement de traction (gauche) et de compression 
(droite) du mortier considere - Comparaison modelisation/experimentation 



Figure 4. Etat de fissuration (mailles de couleur foncee) du mortier numerique sous 
chargement de traction donne par ABAQUS. 


3.3 Effet de la lixiviation sur le comportement mecanique du mortier 

3.3.1. Mecanismes physico-chimiques de la lixiviation 

La lixiviation par I’eau entraine la dissolution de certaines phases solides de la 
pate de ciment hydrate. Le mecanisme a ete longuement etudie et decrit par plusieurs 
auteurs (Adenot, 1999) (Garde, 1996). La degradation est caracterisee par la 
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presence de fronts nets de dissolution-precipitation delimitant ainsi des zones de 
caracteristiques physico-chimiques et mecaniques differentes. La portlandite (CH) 
est la phase la moins stable de la pate cimentaire, elle se dissout en premier. II a ete 
montre par plusieurs auteurs que ce phenomene d’origine physico-chimique entraine 
une chute des proprietes mecaniques : module elastique, resistance a la flexion et 
compression (Garde, 1996) (kamali,1999) (Le Bellego, 2001), angle de frottement 
interne (Heukamp, 2001). Recemment, Kamali et al (2004) ont montre que la 
lixiviation de la portlandite uniquement pent entrainer une chute de module elastique 
d’ environ 50% pour une pate de ciment. Dans cette etude, nous nous interessons au 
comportement global du materiau (zone elastique-pic-post pic). 

3.3.2. Ejfets de la lixiviation sur le comportement en compression et traction du 
mortier 

Seul I’effet de la lixiviation de la portlandite est etudie. La teneur volumique en 
portandite de la pate de ciment est d’ environ 15%. Get hydrate est la phase la moins 
stable de la pate cimentaire, elle se dissout en premier creant une porosite 
supplementaire non-negligeable. Sa dissolution est modelisee en affectant aux 
elements finis la representant les caracteristiques de I’eau. En suivant la meme 
methodologie qu’au paragraphe 3.2, il est possible de presenter les resultats la 
modelisation concernant le comportement en compression du mortier apres 
dissolution de la portlandite. La dissolution de la portlandite entraine une chute 
d’ environ 32% du module elastique (de 42 GPa a 29 GPa) et d’ environ 40% de la 
resistance a la compression (de 63 MPa a 38 MPa) du mortier etudie. Le 
comportement post-pic parait plus ductile apres lixiviation de la portlandite. Ges 
resultats sont en accord avec les observations experimentales trouvees dans la 
litterature (Garde, 1996) (Kamali, 1999) (Le Bellego, 2001). 


Contrainte 

(MPa) 



Mortier sain 
Mortier lixivie 


Figure 5. Comportement en compression du mortier avec et sans portlandite (CH) 
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4. Conclusion et perspectives 


Une approche de modelisation multi-echelles du micro au macro pour la 
prediction et revaluation des proprietes mecaniques et de transfert des materiaux 
cimentaires est presentee. Deux codes de calcul sont utilises, CEMHYD3D pour la 
generation de la micro et la meso-structure et ABAQUS pour le calcul des proprietes 
mecaniques. L’ approche a ete appliquee pour 1’ etude du comportement sous 
chargement de traction et de compression d’un mortier a E/C=0,4. L’effet de la 
lixiviation de la portlandite sur la diminution du module elastique et de la resistance 
a la compression et sur 1’ augmentation de la ductilite sont presentes. Les resultats de 
la modelisation sont en bon accord avec les resultats experimentaux. L’ influence de 
r aureole de transition n’a pas ete prise en compte, elle est une des perspectives de ce 
travail avec revaluation de la diffusivite et de la permeabilite. 
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